Semiclassical theory of current correlations in chaotic dot-superconductor systems. 



P. Samuelsson and M. Biittiker 
Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland. 

We present a semiclassical theory of current correlations in multiterminal chaotic dot- 
superconductor junctions, valid in the absence of the proximity effect in the dot. For a dominating 
coupling of the dot to the normal terminals and a nonperfect dot-superconductor interface, pos- 
itive cross correlations are found between currents in the normal terminals. This demonstrates 
that positive cross correlations can be described within a semiclassical approach. We show that 
the semiclassical approach is equivalent to a quantum mechanical Green's function approach with 
suppressed proximity effect in the dot. 



Over the last decade, there has been an increasing iru 
terest in current correlations in mesoscopic conductors.!!! 
Compared to the conductance, the current correlations 
contain additional information about the transport prop- 
erties such as the effective charge or the statistics of the 
quasiparticles. Systems such as e.g. disordered conduc- 
tors and chaotic quantum dots have been analyzed exten- 
sively, both with quantum mechanical approaches, using 
random matrix theoryaH or Green's functions as well as 
with semifilassical approaches, based on the Boltznum- 
LangevinHQ'Q&El equation or voltage probe models. aau 

Recently, current correlations in normal- 
superconduqting systems have beep-, studied, 
theoreticallya as well as experimentally.!!^ In these 
systems, the current into the superconductor is trans- 
ported, at subgap energies, via Andreev reflection at 
the normal-superconducting interface. In the normal 
conductor, Andreev reflection induces a proxmity effect 
which modifies the transport properties pat energies of 
the order of or below the Thouless energyO However, at 
energies well above the Thouless energy or in the pres- 
ence of a weak magnetic field in the normal conductor, 
the proximity effect is suppressed. 

Nagaev and one of the authorsEl presented a semi- 
classical Boltzman-Langevin approach, in the absence of 
the proximity effect, for current correlations in multi- 
terminal diffusive normal-superconducting junctions, an 
extension of the corresponding approach for purely nor- 
mal systems.Q The interfaces between the normal con- 
ductor and the superconductor was assumed to be per- 
fect. It was shown that the current cross correlations 
are manifestly negative, just as in normal mesoscopic 
conductors In contrast, it was shown very recently, 
taking the proximity.-affect into account, that irL-difFu- 
sive tunnel junctionscfl and chaotic dot juctions,li3 the 
ensemble averaged cross correlations can be positive. 

Interestingly, as shown for the chaotic dot juctionJH 
in the limit of strong coupling of the dot to the normal 
reservoirs and a nonperfect normal-superconducting in- 
terface, positive cross correlations can survive even in the 
abscence of the proximity effect. This suggests that posi- 
tive current correlations could be obtained within a semi- 
classical approach, under the assumption of a suppressed 
proximity effect, if nonperfect normal-superconductor in- 
terfaces are assumed. 



In this paper we present such a semiclassical theory 
for multiterminal chaotic dot-superconductor junctions 
with arbitrary transparencies of the contacts between the 
dot and the normal and superconducting reservoirs. It 
confirmes and extends the result of Ref . [fl5f , providing 
conditions on the contact widths and transparencies for 
obtaining positive correlations. We show that the result 
is in agreement with the circuit theory of Refs. [|l4|l6f| 
with suppressed proximity effect in the dot. This pro- 
vides a simple, semiclassical explanation for the positive 
cross correlations. 




FIG. 1: Upper figure: A schematic picture of the junction. 
Lower figure: Mapping onto an electron-hole junction (see 
text). The arrows show the direction of particle current flows. 

A schematic picture of the system is shown in Fig. [I]. 
We consider a chaotic quantum dot (see Ref. |L7| for def- 
inition) connected to two normal (N\ and N%) and one 
superconducting reservoir (S) via quantum point con- 
tacts. The contacts to the normal and superconducting 
reservoirs have mode independent transparencies T and 
Ts respectively and support N and M transverse modes. 

The conductances of the point contacts are much larger 
than the conductance quanta 2e 2 /h, i.e. NT, MY$ ^> 1, 
so Coulomb blockade effects in the dot can be neglected. 
The two normal reservoirs are held at the potentials 
Vi and V% and the potential of the superconducting 
reservoir is zero. Inelastic scattering in the dot is ne- 
glected. Throughout the paper, we consider the case 
where V\,V2,kT -C A, where A is the superconducting 
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gap. In this case, the probability for Andreev reflection at 
the normal-superconductor interface, R = T|/(2 — T s ) 2 , 
is independent of energy and there-is no single particle 
transport into the superconductor £j 

The semiclassical approach for current correlations in 
normal chaotic dot systems, has been developed within 
a voltage probe approachtltj or equivalently, a minimal 
correlation approach DEEa Here we extend it to a chaotic 
dot-superconductor system, shown in Fig. |l|, by first 
considering the particle current correlations in the dot- 
superconductor system. To visualize the particle current 
flow, it is helpful to view the chaotic dot-superconductor 
system as one "electron dot" and one "hole dot", cou- 
pled via Andreev reflection (see Fig. Q). Due to the 
abscence of the proximity effect, the electron and hole 
dots can be treated as two coupled, independent dots, 
and one can calculate the particle current correlations in 
the same way as in a purely normal system (Andreev re- 
flection conserves particle currents). From the particle 
current correlations, we then obtain the experimentally 
relevant charge current correlations. 

We study the zero-frequency correlator 



P jk = 2 / dt <A/,(t)A/ fe (0)} 



(1) 



between charge currents flowing in the contacts j = 1,2 
to the normal reservoirs, where AIj(t) = Ij(t) — Ij is 
the current fluctuations in lead j (Ij is the time averaged 
current) . The charge current Ij (t) is the difference of the 
electron (e) and hole (h) particle currents, 



I j (t)=If(t)-I*{t), If(t) 



dEif(E,t), (2) 



where the particle currents densities i"(E,t) — if(E) + 
Ai"(E,t). The electron current density flowing into the 
superconductor is i e s (E,t) = is(E,t) = —i s (E,t). 

First, the time average current is studied. Using the 
requirement of current conservation at each energy [due 
to the abscence of inelastic scattering] , the time averaged 
electron and hole distribution functions f e and fh can be 
determined. We have for each dot separately [suppressing 
the energy notation] 



if + *i + ~is = 0, i'l + i$-i s = 0. 



(3) 



The currents i" and is flowing through the point contacts 
are given by 

?? = AT [f a - ff\ , is = MR[f e - f\ (4) 

where ff = [l + cxp([E ±eV j }/kT)}^ 1 is the distribution 
function of quasiparticle type a of the normal reservoir 
j, with — (+) for electrons (holes). Combining Eq. rt^) 
with Eq. (0), we obtain 

fe = MR{ft + /£) + (2AT + MR){ft + fj) 

1 4(AT + MR) ' 1 ' 



and similarily /" = / e ^'\ where we note that f h {E) = 
1 — f e (— E), as demanded by electron-hole symmetry. We 
can then calculate the time averaged currents from Eq. 
(0) and (0). This is however not further discussed here. 

We now turn to the fluctuating part of the current. 
The point contacts emit fluctuations Si" and Sis, di- 
rectly into the reservoir as well as into the dot. As a 
consequence the electron aad hole distribution functions 
aquires a fluctuating part Jail f a {t) = f a + Sf a {t). The 
total fluctuating current in each contact is given by (sup- 
pressing the time notation) 

Ai" = Sif + NTSf a , Ai s = Si s + MR(Sf e - Sf h ). (6) 

The fluctuating parts of the distribution functions, Sf e 
and Sf h , are determined from the demand that the fluc- 
tuating current, in the zero frequency limit considered, 
is conserved in each dot, 

Ail + Ai 2 + Ai s = 0, Ai'l + Ai'> - Ai s = 0. (7) 

Solving these two equations gives Sf e and Sf h , which can 
be substituted back into Eqs. (^) to give the fluctuating 
currents in the leads A?" in terms of the fluctuating cur- 
rents emitted by the point contacts Si" and Sis- The 
charge current fluctuations in contact 1 and 2 is then 
found by subtracting the electron and the hole currents, 



AL = Ai' 



Aij, giving 



Aii = 



(AT + 2MR)[Si\ ~ Si%] + NT(Si^ ~ Si% + 2Si s ) 



2(NT + MR) 



(8) 



and Az2 = A«i^2- For the. second moment, the point 
contacts can be consideredo as independent emitters of 
fluctuations, i.e. (here including Sis) 



{Sq(E,t)Si%(E',t')) = ~8j k S a pS{E - E')S(t - t')Sf (9) 

where the fluctuation power S" (E) is determined by the 
transparency of contact j and the time averaged distri- 
bution functions on each side of the contact £3 as 



Sf(E) 
Ss(E) 



= eAT[/;(l-/;) + / Q (l-/«) 

+ (i-rxr-/;) 2 ], 

= eMR[f e (l - f e ) + f h (l - f h ) 

+ (i - R)(r - f h n 



(10) 



Eqs. (0) to ( |l0| ) form a complete set of equations to 
calculate the current correlations. 

As an example, which also allows a comparison to Ref. 
|l5| , we study the current cross correlations P\2 for V\ — 
V<i = V and kT = 0. In this limit only the energy interval 
< E < eV is of interest, where ff = /f = 1 and = 
f-2 = 0. Correlating the current density fluctuations Aii 
and Ai2 in Eq. dq) and using the expressions for the 
power spectra in Eq. dlfl) , with the distribution functions 
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f e and f h taken from Eq. 



we obtain from Eq. 



and (h ) for the the cross-correlator 



Pl2 = V 



e 3 N 2 MT 2 R 



[2r 2 7V 2 (l - 2R) 



h (AT + MR) 4 
MR(NT[2 -T] + 2MR[1 - T])]. 



fill 



In the limit without barriers, T = 1 and R = 1, as well 
as in the limit r — 1 and N 3> M, this expression co- 
incides with the random matrix theory result available 
in Ref. jl5). From Eq. (pi]), we can make the following 
observations. In the limit of dominating coupling to the 
superconductor, MR 3> AT, the cross correlations are 
manifestly negative (see Fig. |^). The positive correla- 
tions predicted in this limit in Refs. J14|Jl5| ] are thus due 
to the proximity effect. In the opposite limit, dominat- 
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FIG. 2: The cross correlation P12 as a function of the ra- 
tio of the conductances of the point contacts to the normal 
and superconducting reservoirs, 2NF/MR, for different bar- 
rier transparencies. Left: T = 1 and from bottom to top 
R = 1,0.8,0.6,0.4 and 0.2. Right: R = 1 and from top to 
bottom, T = 1,0.8,0.6,0.4 and 0.2. 

ing coupling to the normal reservoirs AT 3> MR, the 
correlations are Pn oc R(l — 2R), positive for R < 1/2. 

In the general case, with different transparencies or 
widths of the two contacts to the normal reservoirs, an 
asymmetric bias V\ ^ V2 or finite temperatures, a de- 
tailed study shows that the conditions i? < 1/2 and a 
dominating coupling of the dot to the normal reservoirs 
are still necessary, but not always sufficient, for positive 
cross correlations. The junction parameters of the exam- 
ple studied above are the most favorable for obtaining 
positive correlations. 

We now show that the result in Eq. (|ll|) can be ob- 
tained by an ensemble averaged, quantum mechanical 
Green's function approach, when supressing the proxim- 
ity effect with a weak magnetic field in the dot. We apply 
the circuit theory of Refs. Jl4|,[l6|], which is formulated 
to treat the full counting statistics of the charge transfer. 
The formulation of the problem is similar to Ref. |Q , so 
we keep the description short. However, we treat arbi- 
trary contact transparencies and finite magnetic field in 
the dot. 

A picture of the circuit is shown in Fig. |^. It con- 
sists of four "nodes", the normal and superconduct- 
ing reservoirs and the dot itself, connected by "resis- 
tances", the point contacts. Each node is represented 
by a 4 x 4 matrix Greens functions (see Fig. ||). The 





FIG. 3: Left: The circuit theory representation of the junction 
with the Greens functions of each node and the counting fields 
shown. Right: The particle counting model, (see text) 



Greens functions of the reservoir nodes are known, given 
by G Nj = e^p(ixjf K /2)G N eyLp{-ixjf K /2) and G s = 
exp(ixfK/2)la x exp(-ixT K /2), where Gat = f z a z + (f x + 
if v )\ and tk — t x & z . Here f{a) are Pauli matrices in 
Keldysh (Nambu) space and Xi > X2 and x are counting 
fields, "counting" the number of electrons entering and 
leaving the reservoirs. The unknown Green's function 
G(xi,X2,X) 01 the dot, normalized as G 2 = 1, is deter- 
mined from a matrix version of the Kirchoffs rule, 



Nl 



I 



N2 



la + h 







(12) 



where the matrix currents Ini,In2 and Is are given by 

W) = NT[G N1(2h G}(4 + T[{G Nm ,G}-2})- 1 
Is = Mr s [G s ,G}(4 + r s [{Gs,G}-2})-\ (13) 

where [A, B]({A, B}) is the (anti)commutator. The last 
current term, dominated by the presence of a magnetic 
field in the dot, isEil Iq, = j(h$/e) 2 [f Gf , G], where $ is 
the magnetic flux in the dot, To = la z and 7 is a constant 
of order of 2 AT, Mrs, he. much larger than one. We are 
interested in the case with suppressed proximity effect, 
which can be achieved by applying a magnetic field cor- 
responding to a flux $ larger than a flux quantum. This 
implies that the 7$ term in Eq. ( |l2] ) is much larger than 
the other terms, and consequently that [f Gfo, G] = to 
leading order in l/$ 2 . The cross correlator Pi 2 is 



Pa = Ve 3 /htr [f K dI N1 /d X 2] 



l Xl=X2=X=0 ' 



(14) 



It is not possible to find an explicit expression for 
G(xi,X2)X) from Eq. ( |i~2|) and the additional condi- 
tions G 2 — 1 and [foGfo,G] = 0. However, to evaluate 
the correlator Pi 2 , only dlNi/dx2 is needed, and con- 
sequently only the Green's functions expanded to first 
order in X 2- We get G = & 0) + i(x 2 /2)G«, where 
& n ' = d n GJdx2\xi=X2=x=0i and similar for the others. 
From Eq. (113) , expanded to first order in X2 as well, we 
then arrive at equations for G^ and &^ . 

The physically relevant result for G^ is f z a z + h(f x + 
i? y )i, where h = NT /(NT + MR). Knowing G<-°\ and 
using that from G 2 = 1 we obtain {&°\G^} — 0, we 
then get &^ = —h 2 ? z a z + h(f x + if y )l + h(f x — if y )l, 



where h = -h[{l -2h + 2h 2 ) + T(h - l) 3 + 2h 3 MR(2R 
1)/(AT)]. Inserting the expressions for G(°) and & 
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into OImi/ d\2, we find Pyi from Eq. fll4|) . This gives 
exactly the semiclassical result in Eq. ([L 1[) 

Finally, it was pointed out that in the limit of domi- 
nating coupling to the normal reservoirs, iVT 3> MR, the 
simple expression for the cross correlations Pyi oc R(l^. 
2R) could be explained by particle counting arguments. E£l 
We now show that the full probability distribution of 
transmitted charges in this limit can be derived from the 
circuit theory. The total cumulant generating function 
F(XijX2tX) is the sum of the functions for each point 
contact, F — F N i + F N2 + F s , where 

F s = eVMTs tr ln[4 + T S ({G, G s } ~ 2)], 
F N1{2 ) = eVNT tr ln[4 + T({G, G N1{2) } - 2)].(15) 

We expand the Green's function of the dot, G, to first 
order in MR/ NT, (but no expansion in the counting 
fields). To zcroth order in MR/ NT, the supercon- 
ducting reservoir is disconnected, and we get directly 
= (drjvi + Gat2)/2. Since normalization implies 
{G (0) ,G (1) } = 0, we note that the contribution to F 
from the point contacts connected to the normal reser- 
voirs, Fjvi + Fm2 oc {&°\ G* 1 '} disappears to first order 
in AT > MR. The total F(xi,X2,x) is then obtained 
by inserting G = into Eq. dTa), giving (for V < 0) 



F = eVM In 



1-R+- (j(xi-x) + e *(xa-x)y 



(16) 



In the limit R <C 1, F has the same form as in Ref. [ p~4| . 
The corresponding probability distribution P(Q,qi,q2) 
of Eq. (|l6|), that gi (92) particles have left the dot through 
contach 1(2) when Q pairs have attempted to enter the 
dot, is (for even q\ +92) given by 



P{Q,Qi,Q2) 



Ql(qi+q2) 



R 



(gi+?2)/2 



(1 - R)Q-(i 1+ i 2 )/ 2 . 



(17) 



<&!©!(Q-[<n+©]/2)!([«i + ga]/2)! 



This is just the distribution function one gets from the 
model of Ref. A filled stream of pairs are incoming 
from the superconductor. Each pair has a probability R 
to enter the dot and then each particle in the pair has an 
independent probability 1/2 to exit through one of the 
contacts 1 or 2 (see Fig. g). 

The processes where both particles exit through the 
same contact contributes oc — R 2 to the cross correlations, 
while the processes where the pair breaks and one par- 
ticle exits throgh each contact contributes oc R(l — R). 
Thus, for R < 1/2, the pair breaking noise dominates 
over the pair partition noise and the cross correlations 
are positive. We emphasize that it is the fact that the su- 
perconductor emits particles in pairs that makes positive 
correlations possible, in a corresponding normal systems 
the electrons tries to enter the dot one by one and the 
cross correlations are manifestly negative. 

In conclusion, we have presented a semiclassical theory 
of current-current correlations in multiterminal chaotic 
dot-superconducting junctions. It is found that the cur- 
rent cross correlations are positive for finite backscatter- 
ing in the dot-superconducting contact and dominating 
coupling of the dot to the normal reservoirs. We have 
also shown that this approach is equivalent to an ensem- 
ble averaged Greens function approach with a suppressed 
proximity effect 
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